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o(h) +0 when h 20), where ¥5 is the solution of the difference boundary 


- ‘ ~h gh yh gh aed 
value problem with disturbed coefficients A,, 3,, De» Fas and u(x) is the °° 
solution of problem (1). Its conservativeness is # necessary and , ee 
sufficient condition for the coefficient stability of the canonical schemes — 
Questioris relating to the convergence and accuracy of the conservative me 2 
difference schemes are studied in §5. It is demonstrated that the : é 
zero-rank conservative scheme converges’ in the class of piecewise . vs 


continuous coefficients (ky a» 29° ); any conservative ‘scheme of the 
first rank has the first order of accuracy in the class of coefficients 


git) (m>1)3 any conservative scheme of the second rarfk which satisfies 
the conditions of the second order of approximation has the second order 


of accuracy in class © e , but in class Q m) (m>1) the first. These 
theorems are proved by means of an a priori estimate W2,f M ipl 9? where 


zis the error in the solution of the difference boundary value problem, | 


and g is the error in approximation of the scheme pike art in the solution 


Cara 7/8 


APPROVED FOR RELEASE: 08/25/2000 CIA-RDP86-00513R001446920010-4" 


"APPROVED FOR RELEASE: 08/25/2000 CIA-RDP86-00513R001446920010-4 


. s/044/62/000/006/075/127 
Uniform difference schemes B168/B112 yee 


of problem (1). By estimating the error of approximation g from the 
norm I, it is possible to reduce the rank of the master functionals 


and the order m of the classes o() or gi@) of the coefficients of. : : 
equation (1). [Abstracter's note: Complete translation.| 7 Ff 
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: $/124 61/000/012/008/038 
a 1d {327 D237/D304 
AUTHORS 3 _ Samarakiyy As Aes and Arsenin, Ve Yao 
TITLE: On numerical solutions of gas-dynamic equations 

with various types of viscosity 
PERIODICAL: Referativnyy ghurnal, Mekhanika, n0Ooc 12, 1961, 

17, abstract 12892 (Zh. vychisl- matem. i 

matem. fiZos 1461, 1, M00 2 357-360) 
TEXT s Pinite-difference schemes of continuous calculation 
(in which lines of discontinuity do not appear) are considered 
for the equations of gas-dynamics of a planes one-dimensional : 
isentropic gas flow with various types of viscosity. Viscosity Vv 


4g subject to the following conditions: The system of equations | 
should have 2 continuous solution, influence of viscosity should 

‘be vanishingly small within the shock layer,and, in the region of 
rarefaction Waves gyugonio ‘Abstracter's note: @ransliteration/ 
conditions should be ‘fulfilled on the boundary of the layer 
puUGONIOT 
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representing a dispersed shock-wave. A rectangular space-time 
net is chosen to formulate difference equations. A difference 
scheme of "cross" type is utilized. The problem of the motion 
of a standing shock-wave is solved, and the solution of differ- x 
ence equations is written interms of moving waves. Problems on 
a standing shock-wave, its dispersion, etc., were calculated 
according to a "cross" scheme, with linear viscosity and Neumann 
viscosity. For linear viscosity, it was found possible to per- 
form calculations with the time jump 3 = 4 times larger than for 
the Neumann viscosity. On the other hand, Neumann viscosity 
guarantees a smaller effective width of the shock layer. Com- 
parison of numerical solutions with exact ones gave fair results 
for the linear viscosity case. / Abstracter’s note: Complete 
translation. 
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of a differential equation of parabolic type 


ian? B168/B112 | oo: “4 
AUTHOR: sanarskiy, ee Ces | si 
aL TLE: A priori estimates for the solution of the difference analog ms | 


PERIODICAL: Referativnyy ghurnal. Matematika, no- 6, 4962, 293 abstract a 
6vi42 (Zh. vychisl- matem. i.matem. fize, Ve 1, 10% 3, 1961, 9 if * ! 
441-460) ia ae 
; ff : 
pextT; In an investigation into the convergence of difference schemes in v 
the class of snooth coefficients, asymptotic orders of approximation and of 
accuracy usually coincide, which is generally not the case with discontinuous fate ad 
coefficients. In this last cases in a study of convergence the application 
of the maximum principle is excluded and the question of estimates of a F 
aifferent nature arises- An a priori estimate in a published work 
(RZbMat, 1960, 40932) is obtained on the assumption of naifferentiability" 
with regard to x of the difference analog of the coefficient of heat 
conductivity. Examining the difference boundary value problem with e 
boundary condi tions of a more general type the author derives an analogous va 
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a priori estimate which is free from this limitation. The a priori 
estimates obtained from the author's integral formula assume that only the . 
grid functions - the coefficients of the equation - are bounded. New i 
a priori estimates are an effective means of proving the convergence and of 
estimating the accuracy of the difference schemes in the class of ; 
discontinuous coefficients. ([Abstracter's note: Complete translation.|:- 
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[6.3909 ! 66500 A060/A126 
AUTHORS : Tikhonov, A. Nes. Samarskiy, A. A» (Moscow) 
TITLE: The Sturm-Liouville finite aifference problem 


PERICDICAL: Zhurnal vychislitel'noy matematiki i matematicheskoy fiziki, ve ts 
no. 5, 1961, 784 - 805 


TEXT : This is a continuation of the work on homogeneous aifference schemes 
reported by the authors in the same journal (v. 1, no. 1, 1961, 5 - 63). The 


solution of the Sturm-Liouville problem for the equation 4 
ploy +Ar(x)u = 0, 04x <1, pleD,y = d_fie(x) G2] ~ a(x)u(x) (1) 


by the method of finite differences has been treated by @ number of authors. They 
treated problems of precision and convergence in the elass of smooth coefficients 
for difference schemes of the partial type. Subject work treats difference schems 
studied in the above contribution, for the nope of the Sturm-Liouville problem 
in the class of discontinuous coefficients Q m), The formation. of the problem, 


the characteristics of the original family of difference schemes, and the proof of 
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the convergence of the finite difference method are presented, With the aid of aa 
a-priori estimate the order of precision for q(m, 1) for the solution of the fiaitte 
difference problem at h-70 is established, It is-proved that the differance 


scheme : 
k,q,ar u 

uf Se My = (ayz), - dy + Apy. 

where fe) is -1 0.5 0.5 
az=[ ( eas! > de J a(x+sh)ds, i = if r(x+tsh)ds, 
-l “0.5 °° “0.5 °° 7 

ensures precision of the second order for the class of discomtimuous coeffisisnis, 
In the continuations to follow the authors promise to treat homogeneous finite dif} 


piece continuous coefficients of equation (1), as well as the problem of precision 
on non-uniform gride, There are 5 references: 2 Soviet-bloc and 3 non-Sovisé- 


bisc, 
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AUTHORS : Samarskly, A, A., Fryazinov, I. Vv. (Moscow) 
TITLE: On the convergence of homogeneous difference schemes for the heat- 


conduction equation with discontinuous coefficients 


PERIODICAL: Zhurnal vychislitel'noy matematiki i matematicheskoy fiziki, v. 1, 
eo no. 5, 1961, 806 - 824 


TEXT: Considerable literature is devoted to difference methods for solving vA 
equations of the parabolic type. A great number of the contributions refers to ; 
equations with constant coefficients. The stability and convergence of schemes — ; 
with continuous, sufficiently smooth coefficients have been studied before by 

others, ‘The case of discontinuous coefficients leads to considerable difficulties 
since, as a rule, in the neighbourhood of a discontimuity the difference scheme 

does not approximate a differential operator. These difficulties may be overcome 

for the heat conduction equation only with the aid of special a-priori estimates. 

In earlier contributions (1959, 1960, 1961) Samarskiy introduced the notion of 
homogeneous difference schemes having one and the same computational algorithm at 

all the points of the difference grid for all the coefficients of the differential 
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equation for some class of functions. In the present work he. considers homoge- 
neous schemes of the skew type for the solution of linear equations of the para- 
polic type with discontinuous coefficients without an explicit distinguishing of 
the lines of discontinuity, 1.e. without any modification of the scheme in the 
neighbourhood of the line of discontinuity for the coefficients. Hence, most at~- xX 
tention is paid to the problem of convergence of skew schemes in the class of dis- 
continuous coefficients, For a quasilinear equation this problem was studied 
earlier by Samarskiy. He proved the convergence of a scheme for the case of mov- 
able ("skew") discontinuities of the heat-conductivity coefficient assuming that 
n°/t 0 as h- 70 andt™ QO. In the present work this requirement is waived for 
the linear heat-conductivity equation, The authors introduce the original family © 
of homogeneous aifference schemes pre and formulate the mixed finite-difference 
preblem,. With the aid of a-priori estimates and the principle of maximum, the 
stability of the schemes as a function of the. initial conditions and of the right- 
hand member is then studied, The case of the. symmetric six-point scheme 


x *# # {0.5 
G * * 0.5 (% = 0.5). 
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is of special interest, For it the stability in the mean value under the initial 


conditions 


and the stability as 4 function of the right-hand member Y: 

lizilo &™ WW i, 
is proved, These estimates hold for the case when the lines of discontinuity of 
the coefficients are the straight lines x = const, and the, coefficient of heat corm 
ductivity k(x,t) and the thermal capacity c(x,t) satisfy the Lipshits condition 
with respect to t. For the case oh a skew discontinuity an a-priori estimate ac~ ae 
cording to the norm of (2° , 1) 1/ , where n = 1, 2, 35 eee The authors conclude 
by giving the convergence proof for schemes of the family under consideration in 
the class of discontinuous coefficients, as well as the estimates for the rate of 
convergence (order of precision) with respect to h and ©, It is proved that the 


scheme ‘f nt with master functionals 
C fe) ¢] 5 


das 4 
aty(s)] = [ ( Hs)! ? pfy(s)] = Fy(s)] = Rty(s)} = {yts)as A) 
-0.5 


“1 
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has greater precision than other schemes of the family under consideration, There 

are 16 references: 11 Soviet-bloc and 5 non-Soviet-bloc. The references to the X 
4 most recent English-language publications read as follows: M, Lees, Approximate 
solutions of parabolic equations. J, Soc, Industr, Appl. Math., 1959, 7, no. 2, 

167 - 1833 M. Lees, Apriori estimates for the solutions of difference approxima- 
tions to parabolic partial differential equations, Duke Math, J., 1960, 27, no, 3, 

297 - 311; M, Lees, Energy inequalities for the solution of differential equations. 
Trans, Amer, Math, Soc., 1960, 94, 58 - 73; H. Keller. The numerical solution of 
parabolic partial differential equations, Math. methods digital computers. N. Y. - 
London, 1960, 135 - 143, 


SUBMITTED; May 14, 1961 
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AUTHOR: —Samarskiy, A. A. (Moscow) 
TITLE : Apriori estimations f 
PERIODICAL: znurnal vychislitel'noy matematiki i 


1, no. 6, 1961, 972-1000 


Ve 


ers difference equa 


TEXT: The author consid 
differential equations o 


a du au du 

ee 2 [ecestou) $e] 1a hm ie i) 
32 

gu = 2[k(e.) 98) — ee t)-ar —! a, t, 


pele.) oe +1(2 


atu 


Gh 
Lau = alk (z, t) =| 


The coefficients k(x, 


functions. 
of the correspondin 


For several cases, 
g difference equa 


f the following (parabol 


t,u) and k(x,t ) are assume 
theerror z=y-U 
tion, is estimated a 
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or difference equations 


matematichesoy fiziki, 


spond . to 
polic) types: 


(1) 
(2) 


tions which corre 
ic and hyper 


250; 


Ou 


t 


=> 
- Ou du 
1 “§23'' Ot 


Ou 
Uy Ge? 


) (3) 
a to be discontinuous 

, where y is the solution 
priori. This 


du 
t,u, Or 
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paper is a continuation of two earlier papers (A. N. Tikhonov, 

A. A. Samarskiy. Odnorodnyye raznostnyye skhemy. Zh. vychesl. matem. i 
matem. fiz., 1961, 1, No. 1, 4-63., A. A. Samarskiy. Apriornyye otsenki 
dlya resheniya raznostnogo analoga differentsial'nogo uravneniya 
parabolicheskogo tipa. Zh. vychisl. matem. i matem. fiz., 1961, 1, No. 3, 
441-460.). A. N. Tikhonov is thanked for revision of the results of this 
paper. Bunyakovskiy is mentioned. There are 13 references: 9 Soviet 
and 4 non-Soviet. The three references to English-language publications 
read as follows: M. Lees. Approximate solutions of parabolic equations. 
J. Soc. Industr. and Appl. Math., 1959, 7, No. 2, 167-183; M. Lees, 
Apriori estimates for the solutions of difference approximations to 
parabolic equations. Duke Math. J., 1960, 27, 297-3113; M. Lees, Energy 
inequalities for the solutions of differential equations. Trans. Amer. 


Math. Soc., 1960, 94, 58-73. 


SUBMITTED: June 24, 1961 
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NIKIFOROV, A.F.; UVAROV, V.B.; LEVITAN, Yu.L.; SAMARSKIY A,, prof., 
otv. red.; ORLOVA, I.A., red.; POPOVA, N.S., tekhn. red. 


[Tables of Racah coefficients} Tablitsy koeffitsientov Raka. 
Moskva, Vychislitel'ndy tsentr AN SSSR, 1962. 319 p. 
(Quantum theory) 
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D299/D303 
1G. 3500 | 
AUTHOR: Samarskiy, AeA (Moscow) 
TITLE: Homogeneous difference schemes for non-linear 


parabolic equations 


PERIODICAL: Zhurnal vychislitel'noy matematiki i matematicheskoy 
fiziki, ve 2, noo 1, 1962, 25 - 56 


TEXT: Homogeneous difference schemes are considered for the non- 
linear parabolic equation 


Pus & (k(x, t) Buy 4 £(x, ty uy gu, Bu) = 0 (1) 


with discontinuous "heat~conductivity coefficient" k(x, t). The XY 
main emphasis is placed on ascertaining the order of accuracy of 


the six-point scheme eae for the third boundary value problem in 


the bounded domain D. Various homogeneous difference schemes are 
simultaneously considered; these are given by means of functionals 
which ensure the second-order of approximation of the scheme, Dif- 
Gard 1/12 
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ference boundary-conditions of the 3rd kind are formulated which 
have -- in the class of solutions of the equation fu = 0 -- the 
same order of approximation as the scheme. The accuracy estimate of 
the obtained difference problem reduces to estimating the solution 
z of a linear difference equation with linear difference boundary- 
conditions and zero initial conditions.» So-called "fixed disconti- 
nuities" of the coefficient k(x, t) are considered, i.e disconti= 
nuities on a finite number of straight lines x = Ny = const., para= 
llel to the t-axis in the (x, +)-plane. It is shown that the third 
boundary-value problem has the same order of accuracy as the first 
boundary-value problem . The main result is formulated as a theorem 
on the uniform convergence of the solution of the difference prob- 
lem. The method used is similar to that cited by M. Lees (Ref. 10: 
Approximate solutions of parabolic equations, J. Soc. Industr. and 
Appl. Matho, 1959, 7, nOc 2, 167 - 183). The results of Ref. 10 
(Op.cit.) follow from Theorem 3 of the present article. The apriori 
estimates used in Ref. 10 (Op-cit.) are however unsuitable for a 
convergence proof with discontinuous coefficients. Par. 1. Homoge= 
neous difference schemes for nonlinear parabolic equations: The 
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problem is formulated as follows: To find a continuous solution 


u(x, t) of Eq. (1) in the rectangular domain D(O<x<1, Otte 
.. Tf), which satisfies the boundary conditions 


k gu - 0, (t)u = u(t), a0 (2) 
x ou ra (t)usuj,(t), x= 1 (3) 
ox 2 2 s im 
and the initial condition = 
u(x, 0) = uy(x). (4) Am 


The coefficients k(x, t) has discontinuities of the first kind on | 
the curves I (V= 1» 25 se+y Vg)+ For convenience, Bq. (1) is re- 
placed by equation 


fu = a (x(x, +) Guy, £0 ty as 2k oa gu = 0, (1') 


which is considered in the following. The problem defined by Ea. 
(1') and related conditions, is called Problem I. The difference 
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equations corresponding to differential equation (1) is obtained by 

means of the integro-differential ‘method. Let OX, = ih, t. = jr, 

i= Q; 15 ooey N, J = QO, 1, eoeg Ly he= V/N, T= T/L) denote the 

mesh, y ~- the mesh function. The sixpoint homogeneous difference 

equation be : . a 

“PE, y= (ays) D + 2.9 (as ty, A yet aY zs Ye) + : (25) ae 
<n + (1— 0) p(x, Ey, ays +642 95) =O, 

is considered, corresponding to Eq. (1'); a (O=<a <1) is a nume- 

rical parameter. The mesh func*ions a@ and ¢ are computed by means 

of the standara functionals 


A[y(s)] (-1<8<0) and F[p(s)] (-0.5<8<0.5), (26) 
vize a= wk (x + sh, t)], 7 = o(x, ts, ¥s avr My + ayzs vz) = 
. F"[£(x Bly hes alty + OY, yz). 
Card 4/12 
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fhe functionals a” ana f° satisfy certain conditions (denoted by D) 
As the original difference system one takes Eq. (25) satisfying con- 
ditions (®). It is shown that the error z = y - u (where y is the 
solution of the difference equation and u —- of the differential 
equation) satisfies a linear difference equation. Thereby the mesh 
function ¥ is obtained, representing the approximation error. Par. 
boundary-value problem: The magnitude of the er- 
irst difference boundary-value problem de- 
The boundary conditions 
d on the difference sche-— 
n the error intro- — 
e difference XN 


My eke’ 3 (u=m(), =o (6) 

[My ke’ + a(t) ual aE (7) 
where u! = du/8x, ti'= 6u/dt). The difference poundary-conditions 
Cara 5/12 
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i) 


are Ay = (al) yx —51y)™ + 0.5h laf (0, t, y, akyxy uz) + (8) 
ee vy 7 a rf 
+(4—a)f (0, 6, ys 2hys, yp) = ul” (t) Iya 2 =0, (8) 
where | = t - te Formulas are derived for the approximation errors 
V4 and Vor of the boundary conditions. The following difference- 


poundary-value : problem is formulated, (corresponding to problem 
(I)): Find the function y = y(x, t), defined on the mesh O and sa- 
tisfying the difference equation 


Fey = (aye) + ap (2, t ys oP ys + avs He) + (26) 
La) plat, ys aie + Stig Yy) =O Sa 9, 
the boundary conditions 


1" dy = ul) (¢), 


1(2)y = ult) (4), 
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and the initial condition = 
y(x, 0) = u(x); (27) 


this problem is called Problem (II). In order to estimate the con- 
vergence and accuracy of the solution y(x, +t) of Problem (II) with - 
respect to the solution = u(x, t) of Problem (I), it is necessa- 
ry to estimate the mesh function z2 = y - u in terms of the ap roxi- 
mation errors , 4, and 5. From the foregoing follows that 2(xXy t) 


is the solution to the following difference boundary-value problem: 
(ag + 0 (2) 083 = — (30) 
z= (az, = 6,2) + hq, (2) — 8:27 =v x=0, (31 ) 
Iz = (azz+ 642)" — hg (2) + $227 = — Va» ee fs (32) 

z(z, 0) =0, (33) 


— 


where Q, q4 and Qo are given by formulas 


Q(z) =.b,42, + Dyo%, + Booty + bo4Zg + dy2 + doz (34) 
Card 7/12 


APPROVED FOR RELEASE: 08/25/2000 CIA-RDP86-00513R001446920010-4" 


"APPROVED FOR RELEASE: 


ECE REET CENEE BRED SARS SON FEA Fae Re ee 


, 08/25/2000 CIA-RDP86-00513R001446920010-4 


33e59 
5/208/62/002/001 /003/016 
Homogeneous difference schemes for --- D299/D303 


qi (2) = (Su2 + Cred 4 Aas2x + Mi2tx) [e=o. \ (35 ) 
qe (z) = (beat 4 Gar + eats + haidz) fete ele 


oem<pcM, - 0<m<acm, | Dx|< Af, In|<M, 
8, > mh > 0, 6,290, 6,+%, 2m >0, [bal <A, 
: |rl<Mo (s, k= 4, 2). 


Problem (30)-(36) is called Problem (III). Further, numerical me- ne 
thods of solution to Problem (II.) are set forth. Par. 3, On the 
accuracy of continued-calculation schemes: Estimates are obtained 

of the solution to Problems (II) and (III). These estimates are de- 
rived in the form of theorems and lemmas. Theorem 1: If the condi- 
tions . 

o 6OM, /o;/ M, o M, /(0,)3/ M, s = 1, 2, 0.5 @ 11,24 


(36) 


s 
are satisfied, then (with sufficiently small t T,) the solution 


to problem (III) is expressed by  __... 
7 "(fest ‘ 


z2(z, t)b< [> th; (Zz, rye] + Tec, HE + fes(z, tih)y< (3) 


=t 
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t=t 
t 


<M {hz (e, 0) b+ jz: (z, 0)b + [2 rp (x, ) lee (3) 


aT 


wnere M and ff are positive constants. On the basis of 2 additional 


conditions 
; Y= ye, YW= we, Va = wy (9) 
and~-~.. : : (10) 
e<M, &<M, &<M ; 


oP SS from the references), one obtains Theorem 2: If conditions — 
44), (2), (9) and (10) are satisfied, one obtains for the solution 

of Problem (III) the apriori estimate wa 
12 (2 <ALL, + FE Db +” (12) er 

vet y v= 12 


[3 ete oni ntel paz wee 


f... ad Pome 

where M = 0 for b,, = Byp = bop = Boy = 0. The solution of Problem 
(III) can be expressed as the sum‘of 2 particular solutions 2 = 2, 
+ Zo, where Z denotes the solution of the problem with homogeneous 
Card 9/12 
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boundary-conditions (v, = Vo = 0), and Zn - the solution with a 


homogeneous equation (W= 0). Lemma 4: For the error in the solu- 
tion to Problem (II), due to the error in approximating the bounda- 
ry conditions, one obtains the estimate 


[eet / awn & t%), 4f 0.5<ae1. (18) 


In the case of continuous coefficients k and f, one obtains AX 


m 
W = o(nh?) + O(t “) at all the points of ©. (20) 


Theorem 3: The solution to Problem (II) converges uniformly to the 
solution u(x, t) of Problem (I), when h and t tend to zero indepen~ 
dently, so that for T < Ty» the inequality 


//y - w//,<Mbn® + 0) for 0.5 <a <1 (21) 


holds if condition (20) is satisfied. Further, the approximation 
error in the neighborhood of the line of discontinuity is estimatc:, 
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With coefficients k and f which are discontinuous, one obtains Theo- 
rem 4: If k and f have discontinuities of the ist kind on.a finite — 
number of straight lines x = Ww = conste, and the (above-) jindica- 
tea conditions are satisfied, then the solution y (x; +4) of problem 
(II) converges uniformly to the solution u(x, *) of problem (1) 
when h and 7% tend independently to 2erd, SO thet for any mesh sys- 
tem 7 ne of the original difference equi bionsy, the estimate 


ify w// 9 < Meal? . 1%) for 0.5<a<1 (50) 


holds; (with sufficiently small tT we) a The author expresses his x 
thanks to AN. Tikhonov. There are 19 references: 14 Soviet—bloc 

and 5 non-Soviet—bloce The 4 most recent references to the English- 
language publications read as follows: J. Douglas, On the numerical 
integration of quasi-linear parabolic aifferential equations. Pacif. 
JeMaths, 1956s 6, nNOol, 35-423 Jo Douglas, The application of sta- 
bility in the numerical golution of quasilinear parabolic aifferen- 
tial equations. Trans. Amer. Math. Soce, 1958, 89, no.2, 484-5185 

M. Rose, On the integration of non-linear parabolic equations py 
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P#RIODICAL: Zhurnal vychislitel'noy matematiki i matematicheskoy sate 
fiziki, v. 2, no. 6, 1962, 1117-1121 18 : 


TEAT: The peopies 


a EX ‘< p : > 
SEA, Lu = 5 (ky (=, t) m=), % NEQr= GrO<t<N, ) 


a=} 
uwlp=us(z, 0), te(0, TI: wie 0) = uo Uo (z), ze6, 
k(t, )2>a>0 (4 =const), a= 4, 2. 


is considered under certain genuine conditions. The domain under 
consideration is covered by a non-uniform rectangular spatial net & 
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AUTHORS ; Tikhonov, A» Ne, Samarakiy, A: Ays (Moscow) 


TITLE: Homogeneous difference schemes with a high degree of 


accuracy over non-uniform nets . 5 


- PBRIODICAL: Zhurnal vychislitel'noy matematiki i matematicheskoy fiziki, 9. 
: v. 3, no. 1, 1963, 99-108 _s 


TEXT: This paper is a continuation of previous: papers!’ (Zh. wehisl meted: YA : 
matem. fiz., 1961, 1, No 1, 5-64; and No. 3, 425-440), the fundamental. °/% ‘ 


estimates of which are shown to be valid without the additional conditions: a 
C _ 
OCMC Hy /HyK My» AYR Boe. (2) 


Phe accuracy of zero-rank schemes is characterized by the mean square 


step 


; ; N 
tine = Cn2)'f? ow D remy! 
He jot ivi . 


It we to o( fill 3). The results of another paper by the same authors 
Card 1/2 : a } : 


APPROVED FOR RELEASE: 08/25/2000 CIA-RDP86-00513R001446920010-4" 


"APPROVED FOR RELEASE: ei 25/2000 CIA-RDP86-00513R001446920010-4 


y 


s/208/63/005/001/004/013 
Homogeneous difference schemes with ... B112/B102 


(Zh. vychisl. matem. i matem. fiz., 1962, 2, No. 5, 812- 832) concerning. 
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: <. pam awe ere ee ee ee Ba i ee oe uae, 
(TLE: ss: Schemes of uniform differences over uneven lattices for parabolic~_ 
type equations a oe Me ae ee eae eae 
.» PERTODICAL: Zhurnal vychislitel'noy matematikt 1 matematicheskoy fiziki, Vo ee 
7 nO. 2, 1963, 266-298 (9 es eee ae Bes 
TEXT: | tn previous papers (Refs 1-5: Zhe yychiel. naten. i matem. fine, 1961) 
1, Now 3, 441-460; 1961, 1, No. 5, 806-8245 1961, 1 Noe 65. 979-1000; 1962, 2, Now 
1, 25-565 1962, 2, Now 45 603-634) .the author investigated ‘the basic problems. of 3 


- ‘the theory of uniform differences of ‘Jattices for linear, quasi-Linear, and non=" 
m linear parabolic-type equations. and established the. stability; convergence, and. - 


_ estimates of the rate of convergence of certain families of uniform “differences -< 
within the class of continuous and discontinuous: coefficients of differential equa 


 phorige fe Ne Tikhonov and-A. A. Samarsiiy (Ref. 6: Zhe vychisle matem. i matems” 


” Bize, 1962, 2, Now 5; 812-832) pointed out that schemes of differences having & 
second order of approximation over an even lattice have only & first “order of ap- 


_. proximation over uneven lattices indicating that the acciracy problem for uneven 
. dattices requires further studies. ‘The present paper. studies schemes of uniforn . eS 
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ptksas2), 2 z (x(x) $ 2) “a(x)a + @) = = oO 


which over special k, q, and f-dependent sequences. of lattices a (8) have ‘for: 


. the class of discontimous. coefficients a second degree ‘of accuracy... “The studies i 
-of the errors, discontinuous ‘coefficients, - and order of. magnitudes over uneven -: a 
- lattices are followed by the investigation o of uniform Schemes, of " differences for 
. the linear and quasi-~linear equations , : eee 


e(x,t) ge - 2 Poe (e(x,t) 2 rae + (est) 3 ae =. -aCxst) u+ + Flat) pon 


otest) 3 an 2 (bast, ak + #0, 1 ‘ty. y ge e | 2 Se. om 


_ and the author caleulates their: respective. errors ‘of approximation oo sr : a 
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priori estimates and uses then for the ere = 
On a Sequence of uneven lattices, 
these schemes remains unchanged during the trandition ke the order of ‘accuracy of 
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| TITLE: Locally one-dimensional difference ‘achenes on nonuniform Beta, Pe 


SOURCE: Zhurnal vychtelitel !noy matematikt i natonatichestoy fiat, Va: aan eo 
no. 3, 1963, 431- 466 — Rote ; 


TOPIC TAGS : parabolic-type differential equation, dtttecence scien ae 


“ 


ABSTRACT : Difference gcheuier ay the solution of: _prdimenatonal parabolie-type 


a ee 


equations Bae ee ; ; 
olx,)2E = : Ss Let +f, o2 fae Bs 
oes ee - Ged, oe a Te 


: where x = Gegetenie ” is che potae’ of ‘p-dimensional Buclidean’ ‘space, ° L. ts ‘a 
‘differential ‘operator, and c(x,t) and £ are. certain functions, are studies for: 
~ four particular ee of ne satekwscpes operat or lg ans the. function’ £ o-: 


oo cole nena anne Tiger eet meedntinn Caine 
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ACCESSION NR: AP3001099- 


. Depending on the forn. of Lag and £, fae secu léus ly (k.= 0, a: 3) are’ - defined 
for the solution of which the construction of. ‘locally. one-dimensional difference: 
schemes IIk (k = 1, 2,:3, 1!, 21) using nonuniform nets is presented. nee 
Determination of the approximat tori: error z(x, t)-y(x,t), where u(x,t) is- the. 
solution of ‘the problem I, and y(x,t) is the solution of the corresponding | 
difference problem TI, characterizing the correctness: of the schemes, -is : 
analyzed. Great attention is paid to the a’priori estimates (mean. estimates and: 
uniform estimates) of the correctness of schemes iI,. To obtain such estimates 
the method of energetic inequalities of nth rank developed earlier by the author 
is used. The a priori estimates derived make it possible te prove the — 
convergence of locally one-dimensional schemes when the coefficients.of the — 
differential equation have discontinuities of the first kind on a Finite number | 
of hyperplanes parallel to coordinate hyperplanes. It is noted that by selecting : 
proper nets, second-order accuracy with respect to h.(h is a mesh size) canbe = = 
achieved for the difference schemes -Paeannreds (Orig. arte eat 168 waver teats ear) 


“ASSOCIATION: none. 9 0 pans 
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‘AUTHORS: Samerskly, As Ao Sobol', I. NM, (Hoscon) a aa i : . | 
“TELE: _ Exemples of numerical computation of te rature wave), \ 

as coupora\ure Wave 

t 


; : oe Zhurnal vytchisl. matematiid 1 matonattote fia, v. 3s no. hy 298, 702: bas. 
(9. 


“ In0Ptc TAGS: differential 1 equation, | ners equation, generalized eer ee say : 
"(solution 2 ee 
no ' ABSTRACT: This paper 1s concerned with numerical solution of a cuast linear cine ; 
a ! of, heat conductivity _ cata ite Satta ’ fade Ber 


‘for the cases Pp = (4,2)3- b ‘usual it is always assumed that : Ee 
im 
| 


Ka(u) «= xu", 


| 
pe ale 
“ cs cee 
‘where oy zl, > 0. 0. athough (Q) arises in various areas: of mathematical physio, | 
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ad pee 


-- lene authors, for definiteness, call the function u = (t,x, 0+0s%p) a temperature i 
function. Ya. B. Zel'dovich and A. S. Kompaneyets (K teorli rasprostreneniya-tepl# | : 
lord teploprovodnosti , zevisyashchey ot temperatury*. Ngh..k semidesyatiletiyu - 
‘lakademika A. F. Toffe". He, Izd-vo AN SSSR, 1950, @i-71.) and G. I. Barenblatt (0 |. 
~ :inekotory*kh neustanovivshikhsye dvigheniyakh zhidkosti i gaze v poristoy srede. - 
- Prikl. matem. i mekhan., 1952, 16, Ko. 1, 67-78.) have shorn that equation (1), in |. 
‘case -p = 1, has a solution whose derivatives, et the points where u(t,x) goes to} - 
igero, are discontinuous and the flow K(u) 2 u/9 x is continuous, 4.0. there exists | - 
. ta temperature front u = 0 which is propagated with finite velocity. 
solution. The existence of a generalized solution of 
ue problems are proven by 0. Ae Oleynik, A. S-: 
pa nestatsionarnoy filttratsii. Izv. B.~ nee 
oved convergence of an explicit ye 


H 
“t 


a) ae tee 


can probably be extended ‘to ‘the é 


‘in the class of generalize 
- ease of imolicit schemes) odnoy zadachi nestatsionerroy . ib 


filttretsii metodom integr vytchisl. natem. 4 matem, fizes 
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11962, 2, No. 5, 938-942) computed the generelized solution of en equation of the - |: + 
‘form (3) by Dorodnitsy*n's method of integral relations. For computation of svch: , 
‘generalized solutions (which. they call temperature waves or simply solutions) the: 9! .! 

" “tanthors use homogeneous difference schemes of continuous computation not specifying.! . 

' ‘elear separetion of points of weak discontinuity. The theory of such schemes has“: . 

_- been worked out by various authors. However, all the proofs of convergence assume Pe 

. that Ky (vu) 2e>0 and despite the great generality of these theorems, they sre not 


' ‘applicable to the case where Ky (u) goes to zero (even allowing discontinuous - 


. _ functions R= Ky (t,x,u)). ‘The aim of this article is to show that these schemes! — : 


‘are also suitable for computation of temperature waves. Such schemes make it: 9 i». 
‘possible to carry out the computation by large steps in time, to give the velocity jo: 
of propegation of the front well, and, for a sufficiently fine grid, also the” : 
‘profile of the front. In the case of space variables (py) the authors use the .* 
‘locally one-dimensional method of variable directions set forth by A. A. Samerskiy { . 
- '@% ognon ekonomictmom raznostnom metode reshoniya mogomernogo parabolicheskogo 
- ~wravneniya v proizvol'noy oblasti. Zh. vy*chisl. matem, i motem. i fiz., 
| 9062, 2, No. 5, 787-811) and (Local'no-odnomerny*ye skhemy* na neravnomernoy bo 
“getke dlya mnogomerny*kh psrabolicheski kh vravneniy. Zh. vy*chisl. matem. Zi metem. j7! 
fine, 1963, 3, Ko. 3, 431-466). The authors give a brief charecterization of the | | 
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. |ACCESSTG! NR: AP3004958 a5 Se a4 te 
‘method of the last two cited pepers which is applicable to equetion (1). A step in, 
- letme +) 2 t—< ttl is divided into p layers of identical thickness ( "fractional 


x 4 steps") 


> oe 


nee eee 


| ; ee ae eee i | 


‘ltn the leyer numbered 2 one solves the one-dimensional equetion po a 


I a ee 

: | | ‘ 

, -|Here, all other coordinates Xp» distinct from x, » play'the role of parameters. |}. 

|. tAbthis stege, for boundary conditions t ¢ the boundary func- ‘| 


(aes at points of intersection of stra axis.from’- 1°. i. 


a 


mogeneous difference sc 
pinion: there is no more sui 
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‘the solution of ear parabolic equations. This method is 
- lapplicable to arbitrary reg ( to parallepipeds) and keeps its order of = 
 laccuracy on inhomogeneous grids.. 4 linear ‘parabolic equetions > 
of general form even in the presenc continuities (type I). 
‘In..such a wide area of usefulnes 3 a whole 
be jsories of merits; simplicity of 
- jof other schemes) of the requi 
' : tion with very coarse steps in time. 
_irapidly complex problems where great accura 
: ‘difference scheme give, instead of the exac 
_l(the finer the grid, the greater the accurecy) - 
_ jthis profile with a very coarse grid end estima 
«ithe euthors constructed, in§5, for the case p 
-tanalog to the well known equstion of the form u 
‘\(the constant ¢ is the velocity of the wave). °F 
‘‘conputation of gas dynemics with viscosity the difference running wave was con- + 
-igtructed by A. A. samarskiy and V. Ya. Arsenin (0 chislennom reshenii vramneniy. 9. 7 
‘pazodinamiki s razlichny*m tipemi vyazkosti. Zh. vy*chisl. matem. 1 metem. fiz., Hee ape 
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ABSTRACT: The study of traveling waves under conditions of nonlinear theraal con 
ductivity is related to the problem of a piston operating under thermal conditions. 

In the framework of the one-dimensional ptane problem for hydrodynamic equations , 
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= i dance of this width on the parameters of ‘the problem is established. — 
’ | oF the traveling wave is regarded as a model for the analysis of poss t 
‘. i and their dependence on the degree of nonlinearity of the heat equation. The field 
lof the integral curves is analyzed in relation to the values of indices of non~- 
| linearity. In a number of cases, analytic solutions are presented. A classifica- 
‘ticn of various types of traveling temperature waves is given. Difference sethods 
are used for the machine solution of a system of partial differential equations 
for apt ropriate boundary conditions on the piston, and the results for a number of 
_ : computer solutions are presented. A comparison of the analytic results with nu- 
* | merical solutions makes it possible to judge the accuracy of the difference methods 


‘used and to affirm the existence and stability of the traveling waves constructed. 


ea 
- |"The authors are grateful to L. N. Busurina and V. P, Krus for programming and per~ 


| forming the computer calculations, and also to Loy Lok'yanova, A. 4. Zakharova 
“and N. G. Teplova, who participated in the work of doing the calculations end 
formulating the results. The authors acknowledge the useful discussions offered by 


B. L. Rozhdestvenskiy." Orig. art. has: 27 formuias, 14 figures. 
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‘AUTOR: * gamer akiy, A, A. (Moscow) g 2 


TITLE: On monotonic ‘difference representations for elliptic and perabolioc equations! 
jin the case of a non-gelf-conjugate elliptio operator 


SOURCE: Zhurnal vychislitel'noy matematiki i matematicheskoy fisiki, v. 5, no. 3, 
1965, 548-551 


va TOPIC PAGS: Laplace tranaformation, “elliptio: equation, parabolic equa ag ti on : ee 
| [@ppracimation method vA 


ABSTRACT: Monotonic representations for a non-self-canjugate elliptic equation of 
jsecomi order in an arbitrary region are developed. The first schemes considered are 
|charaoteriged by a uniform rate of decrease denoted by O(h2). Additional considera-' 
ition is given to monotonic, locally one-dimensional schemes shich decrease uniformly 
‘at a rate O(h2 + T ) for a parabolic equation with a ron-self-conjugate elliptic _— 
‘operator. The monotenio representations are demonstre.ted by means of example 
‘problems. The first example yields . representation cf the firat boundary problem . 
in a particular region of a multidimensional space. The problem ia defined in - F 


terms of ite first derivatives, end a monotonic representation is developed to the 
Card {| . 
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= secon order: of accuracy. she ‘author’ proves the validity of the presumed: eseareay: 
of the scheme developed. The second example involwes a parabolio equation in the 
aylinirical region r= (@+1) xKo<tch. fhe problem ia stated as 
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Additilonal conditions and definitions are 
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The selection of a time Sitervel is discussed, and certain substitutions are unde 
auch that the locally one-dimensional acheme is represented in the form 
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Certain ramifications of the problem due to partioular boundary oonditions | 
| discussed. Orig. art, has: 11 equations. . ! 
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TITLE: Theory of difference schemes 
SOURCE: AN SSSR, Doklady, Ve. 165, N0. 55. 1965, 1007-1010 - 
TOPIC TAGS: Cauchy Froblen, Hilbert epace, Banach space 


ABSTRACT: _ The abstract Oauchy problem ee a ea ea ee 
da /dt + A(t) = 1(2), Ostet, u(0) = = Uo, moe D(A) “() 


New ee 
in Banach and Hilbert é spaces is éonsidered.. “The neh permits. determination of. 
sufficient conditions for correctness of ‘schemes which make use only of operator = | = 
properties which are simple to verify directly, These include such properties as 
semi-boundedness from below, etc. The author expresses his gratitude to A, N. 
Tikhonov for his valuable discussions of the results, This paper was presented y 
academician M. V. Keldysh on 19 April 1965, Orig. art. has: 10 formulas. — 
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TITLE: An economic scheme of through ‘Galotius for imultidimensional Stefan © 
probleme te Re Coe eae ae ae 


! SOURCE: Zhurnal vychislitel 'noy matematild 1 -matematicheskoy figiki, ve 5; nos 5, 
; TOPIC TAGS: differential equation, heat tranefer, calculus 


| ABSTRACT: The authors propose an economic difference scheme of through calculus <7). 
' for the numerical solution of Stefan problems with several spatial variables and =: 
several numbers of phases. In the proposed scheme the boundary of phass division 
was not chosen explicitly and homogenous difference. schemes were utilized. The 
proposed scheme was tested for the case of one and two spatial dimensions. In all! =: 
cases the method yielded solutions with adequate accuracy, The authors solve the... 
multidimensional quasi-linear equation of heat transfer by employing the locally =.” 
one-dimensional method. This consists of e step-by-step solution with respect to 
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various spatial variables of the ogsdtnend tiene: heat transfer equation ‘via ae be | 
unconditionally stable implicit schemes. ‘The authors express their thanks to. i 
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I. Ve Friazinov for advice. Orig. art. has: 3 figures, 20 ea and 6 tables.!: 
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TOPIC TAGS$ algorithm, differential equation, computer technology, mathematic phys: ¢ 
ABSTRACT: The article describes difference scheme s in somneeh ar with the development 
of algorithms of a universal t for computers, suitable for sate fares) 
of problems /with the same programs. Orig. art. has: 66 formulas. 
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POPIC TAGS: Cauchy problem, approximate numerical 
solution, mathematical model © . aoe 
olution of multip! e-dimensioned 8 : 


ABSTRACT: Locally one-dimg signal systcus of § 
problems of mathematical shysivs are studied. The method presented: is pased upon. | 
the principle of modeling & muti-dimensioned differential equation. with. the aid of. 
a system of one-dimensional 4q0htions with successive application of simplified oes 
difference systems to the solution of each of the one-dimensional equations. The 
‘principle is applied by the author to the abstract Cauchy problem expressed aS 
du/dt + A(t)u=/(2), o<ixl, u(0) =!t0, nao 
‘ eH and 


‘where A(t). is a iinea® unbounded operator in space H, and u = ult 
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f = f(t) & H are defined on the interval O<+t<¢. Of interest isthe case where - 
A is the sum of linear operators Key ; ae eon i : 


! “A= >» Ae ©. ) 

' A total of m one-dimensional problems: is ‘ised? in modeling the: milti-dimensional 
problem. .In each interval {ts, ti+1] the system of differential equations 
dyya)/ dt + Aa(t)%2)(t) == fa(t), Re 3 
te (4tH], 7 =0,4,. oe 

vis successively solved with initial conditions” 

Yay (#8) = Ye-y(t4), a > 4; oS om (t8) == v(H);: 

720, v(0)= 

Approximate solutions to this system are ee enaae for varying conditions of ; 
parameter bounds. Five theorems are developed for several cases, and additonal — 


cases are referred to earlier research. This paper was presented by Academician 
M. V. Keldysh on 19 April 1965. Orig. art. has: 8 equations. , 
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TITLE: Some problems in the theory of difference schemes 


SOURCE: Zhurnal vychislitel'noy matematiki i matematicheskoy fiziki, v. 6, no. 4, i 
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WW. 
TOPIC TAGS: numerical analysis, difference sobwme, 1 “7 3 ’ ; 
scheme=stebility,__scheme-eonrergance : oa a roll ne shes egpntio : 
ABSTRACT: The algebraic identity and stability of longitudinal-transversal schemes ; 


(special form of locally uniform sckemes) are analyzed by means of iterative schemes 
for solving the stationary equation 


(AtAyv=h? 


where A, and A, are linear positive-definite operators in a real Hilbert space, It 
is shown that longitudinal-transversal schemes are algebraically identical to three 
simple,locally uniform schemes, and new estimates of the rate of convergence of such 
schemes are derived. It is shown that in order to decrease the initial error i/e 
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times (e > 0 is the foreassigned accuracy of the scheme), using the longitudinal- 
transversal scheme, it is necessary to carry out 


“y@)>e yl | (2) 


iterations (n = 6/A, where 6. and A are certain positive constants). The derived 
estimates are applied to the solution of three problems: 1) the difference Dirichlet 
problem in an arbitrary domain for the self-adjoined elliptic operator with variable q 
coefficients; 2) the Dirichlet problem for the nonself-adjoined elliptic operator : 
with constant coefficients; 3) the integro-differential equation of radiation trans- e 
fer. The longitudinal-transversal scheme for solving the difference Dirichlet 
problem in an arbitrary p-dimensional domain applies in the case where the difference 
Laplace operator is represented as the sum of nonself-adjoined operators. Orig. art. 
has: . 90 formulas. , , {LK] 
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92-58-3-15/32 
AUTHORS : Anisimov, A.F., Engineer and Samarskiy, A.G., Engineer 
ee 
TITLE: Continuous Control of Neutralization of Crude Stock 


in Atmospheric-Vacuum Pipe Stills (Nepreryvnyy kontrol' 
rezhima podshchelachivaniya syr'ya na ustanovkakh AVT) 


* PERIODICAL: Neftyanik, 1958, Nr 3, pp 15-16 (USSR) 


ABSTRACT: To protect the equipment against corrosion the crude 
stock processed in an atmospheric vacuum pipe still nas 
to be neutralized. The neutralization of crude stock is 
ordinarily controlled by a laboratory test of the 
water settled from the wide-fraction in the surge tank. 
The lack of continuous control over the alkalinity of the 
water shortens the service life of the processing 
equipment. Therefore, the special construction bureau 
for the introduction of automation in refineries has 

studied this problem at the atmospheric-vacuum pipe 

still of the new Ufa Refinery and ascertained that the 
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pH value of the condensation water settling in the tank of 

the evaporator is almost always 10, while the pH of water 
drained off from the surge tank of the atmospheric tower 
depends on the concentration of the reagent (Fig. 1). As a 
result of this study, it has been proposed that the pH value 

be determined in the stream before the stream enters the surge 
tank. Therefore, samples of the product mixed with the settled 
water were taken at the lowest point of the stream line section 
between the condenser-cooler of the atmospheric tower and the 
surge tank (Fig. 2). The mixture was directed through a 1/2-in. 
tube to a small settler, where the separation of water from 
the mixture took place. The settled water was filtered and 

its pH determined. In the course of this operation, it has been 
found that the present system of calomel electrodes is inadequate 
for determining the pH value in a polluted mixture containing 
petroleum products. Therefore, the calomel electrode design was 
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modified (Fig. 3). As 2 result, the redesigned electrode can 
now be used during « much Longer period of time, and can be 
easily dismantled and reassembled, On the whole, the equipment 
operates satisfactorily and this fact wag confirmed in 1957 by 
the new Ufa Refinery. There is one diagram showing alkali con- 
sumption, one drawing showing the continuces determination and 
recording of the pH value, and one sketch of the redesigned 
calomel electrode. 


ASSOCIATION: Ufimskoye ctdeleniye konstruktorskogo byuro po 
avtometizetsil (Ufe Section of che Design Bureau 


for Automation) 


AVAILABLE: Library of Congress 
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$0: Letopis! Zhurnal'aykh Statey, No. 29, Hoskva, 1949 
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"Vertical Gradiometer," by M. V. Samarskiy, Nauk zap. Luts'k. 
ped. in-tu, 3, No 2, 1955, pp 1Ol-114 (Ukrainian) (from Refer- 
ativnyy Zhurnal -- Astronomiya Geodeziya, No 11, Nov 56, Ab- 
stract No 6554, by G. S. Razdymakha 


"The possibility of determining, by means of an instrument called 
a vertical gradiometer, of the second derivative W of the potential 
of the force of gravity with an accuracy of 10 etvésh sufficient for the 
survey and compilation of anomalous masses located under the observation 
points at small depths and producing an effect of several tens of etvesh 
is presented. For this end we may use a gravimeter with a response up to 
0.001 mgal, equipped to record two values of gravity accelerations at the 
same observation point with two installations of the gravimeter in two po- 
sitions with a constant elevation difference of one meter. During the 
short time period necessary for observation the error of computation of _ 
zero-point variation of the gravimeter, the temperature variation in the - 
thermostat and the barometric correction will be eliminated. 
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